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Abstract 

We demonstrate that by making satellite-based parallax measurements of Ma- 
cho events, it is possible to distinguish disk from bulge Machos and to determine 
individual masses of disk Machos to an accuracy of 0.2 in the log. This is at least 
a 3-fold improvement over what can be done without parallaxes, i.e., just from the 
observed time scale of the events. In addition, we show that the physical distribu- 
tion of disk Machos can be found from 'reduced' velocities determined by parallax 
measurements. 
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1. Introduction 



Major efforts to detect dark matter by observing gravitationally lensed stars 
located in the Galactic bulge are being carried out by the OGLE (Udalski et al. 
1994) and MACHO (Alcock et al. 1994). Several dozen candidate events have 
already been reported. The light curve of a lensing event is given by 

M«t)]= " 2 ^ 1/2 , «(*) = vW" 2 ('-'o) 2 , (i.i) 



+ 4)V2 : 



where u; -1 is the characteristic time, (3 is the dimensionless impact parameter, and 
to is the midpoint of the event. Of the three parameters in equation (1.1), only lu^ 1 
provides information about the lenses responsible for the events. The characteristic 
time is related to the physical parameters of the lens by 

_! (4GMD ol D ls /D os f 2 
U ~ Vc ' (L2) 

where M and v are the mass and transverse speed of the lens, and D i, D os , and _D/ S 
are the distances between the observer, lens, and source. However, these physical 
parameters cannot be measured separately, and hence the physical information 
about the lenses is poorly determined. 

There have been two general approaches advanced for finding additional infor- 
mation about individual events. The first one is to measure space-based parallaxes, 
which yield "reduced" transverse velocities v and corresponding "reduced" Einstein 
ring radius R e (Gould 1994b, 1994d), where 



and 

When a single satellite is used for the measurement, there would generally ex- 
ist a four-fold degeneracy of v; two-fold by magnitude and the other two-fold by 
direction. This degeneracy can be completely broken by the measurement of par- 
allaxes from a second satellite (Gould 1994b). However, for most Macho events a 
single-satellite parallax is sufficient to measure v and R e , leaving only a two- fold 
degeneracy in direction (Gould 1994d). The other approach is to measure the lens 
proper motion fl = v/D i by photometric (Gould 1994a; Nemiroff & Wickramas- 
inghe 1994) or spectroscopic (Maoz & Gould 1994) methods. 

In this paper we demonstrate that the events generated by bulge and disk 
Machos could be distinguished and the individual masses of disk Machos could be 
measured with a significantly improved accuracy from the reduced speed measure- 
ments of events seen toward the galactic bulge. Since the required parameter is the 
reduced speed not the velocity, this method eliminates the need to launch a costly 
second satellite for the complete measurements of v. We also discuss a method 
to extract information about the physical distribution of disk Machos from the 
measured reduced speeds. By contrast, we find that it is only rarely possible to 
measure proper motions of Machos seen toward the bulge. Moreover, even when 
they can be measured, proper motions in this direction provide less information 
about the Machos than parallaxes. 

The reason that parallax measurements are sensitive to the mass of disk Machos 
while proper motions are not, can be understood as follows. The flat rotation of 
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the disk implies that the mean transverse velocity as a function of distance, v(D i), 
is given by 

v(A^) = §^o,0), (1.5) 

where vq = 220 km s" 1 is the rotation speed of the Sun and Ro = 8 kpc is the 
solar galactocentric distance. Hence, the mean reduced speed is a definite function 
of D [. Measurement of v therefore strongly constrains D \. To the extent that 
v and D i are known, the mass is also known. However, from equation (1.5), the 
mean angular velocity Q(D i) = (vq/Ro,0) is independent of distance. Hence, a 
proper motion measurement provides almost no information about the distance. 



2. Probability Of Lensing Events 

2.1. Bulge Lensing Event 

The optical depth to lensng is defined as the angular area covered by the 
Einstein ring A = n(R e / D i) 2 as a fraction of the sky. If the bulge stars are 
distributed over distances between d\ and di from the Earth with a number density 
distribution n(D 0S ), the optical depth r of a bulge star lensed by foreground bulge 
stars is 



4ttC7 

Tbulge — ^2 



dD os n(D os ) J dD ol p(D l)D 



dD os n(D os ) 



(2.1.1) 



where D = D [Di s /D 0S , and p(D {) is the mass density of Machos along the line 
of sight. As the distance from the observer increases, the volume element and 
total number of stars in the element increase. However, the actual number of stars 



one can detect decreases with increasing distance due to detection limits. In our 
computation, we assume that these two factors cancel out. This is equivalent to 
the P — — 1 model of Kiraga & Paczyhski (1994, hereafter KP). 

The frequency of lensing events lensed by the bulge stars, r& u / 5e , is 

D i 



r bu i ge = ^^ J dD 0S n(D 0S ) J dD olP {D ol )D 1 / 2 



J dv y J dv z vf(v y ,v z 



d 2 - 1 



dD 0S n(D 0S ) 



(2.1.2) 



Ldi 

where v y and v z are the components of v and f(v y ,v z ) is their distribution. Here 
the coordinates (x, y, z) have their center at the center of the Galaxy and x and z 
axes point to the Earth and to the north galactic pole, respectively. 

2.2. Disk Lensing Events 

The source bulge stars are located within a narrow region compared to the 
typical values of D i and D[ s of disk lenses. Therefore, it is possible to approximate 
the source stars as being located at a fixed distance. We assume D os = 8 kpc. Then 
the optical depth is 

d 2 

Tdisk = ^rj dD ol p{D ol )D. (2.2.1) 
o 

Similarily, the event rate of a bulge star being lensed by the disk stars, Tdisk-, 

is 

Dos 



^disk=^^ J dD olP {D ol )D 1 / 2 J dv y J dv z vf(v y ,v z ). (2.2.2) 
o 
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3. Models 



To compute r and T, we adopt various models of velocity and density dis- 
tributions. The results are used to investigate the sensitivity of r and T to the 
choice of model. We always assume that the number density of the bulge sources 
is proportional to mass density of the bulge lenses, nap. 

The transverse velocity v is defined as 



V/ - Vlos = V/ - 



Dai , D ls 
v s — + v — 



(3.1) 



where v/,v s , and v are the transverse velocities of the lens, the source, and 
the observer, respectively. We assume that the velocity distribution is gaussian 
f(v y ,v z ) = f(v y )f(v z ), where 



f( v v) = i ex P 

fa** 



(3.2) 



and similarily for f(v z ). We adopt v Zjdi8k = v zMge = and a Zjdisk = 20 km s \ <T Z:bu i ge 
100 km s _1 for the z component of velocity. For the y direction, the values 
adopted are v y ^i s k — 220 km s" 1 and <J y ^ isk = 30 km s^ 1 for the disk lenses, 
whereas Vy^ulge = and cr yj bulge = 100 km s _1 for the bulge lenses. Another 
disk model where the velocity dispersion increases steadily toward the galactic 
center (Lewis & Freeman 1989) is separately considered. In this model (lin- 
ear disk velocity model), the velocity dispersion in the solar neighborhood is 
(°y, disk i disk) = (30,20) km s _1 , and it increases linearly toward the galactic 
center until it becomes y , disk ■,< J Z , disk) = (75,50) km s _1 at the galactic center. 



The influence of bulge rotation is tested separately by assuming that 



Vrot — V-max 




(R < 1 kpc, rigid body rotation), 



Vrot — "Umax 



- (R > 1 kpc, flat 
R' 




(3.3) 



where v. 



max — 



100 km s 1 , and R = \J x 2 + y 2 . 



We also consider the possibility of a bar-like structure for the bulge and its 
influence on lensing events. In the model the bar has axis ratio (1 : 0.33 : 0.23) 
and it is rotated in the galactic plane with its near side in the first galactic quad- 
rant making an angle 9 = 20° between its major axis and the line-of-sight to the 
galactic center (Dwek et al. 1994). The velocity distribution for the barred bulge 
is deduced from the tensor virial theorem (Binney & Tremaine 1987) with result 
(oy, ay, oy) = (113.6,77.4,66.3) km s" 1 . Here, the coordinates (x\y',z') have 
their center at the galactic center, x' represents the longest axis, and the shortest 
axis (z ; ) is directed toward the north galactic pole. Because the mass-to-light ratio 
of the bulge is not known a priori, the virial theorem gives only the axis ratios of 
the velocity dispersion tensor. To obtain our result, we normalize to the observed 
mean line-of-sight velocity dispersion of ~ 110 km s _1 . This calculation yields a 
total bulge mass Mt> u ig e ~ 1.2 x 10 10 M Q , very close to Dwek et al.'s estimate of 
Mbuige ~ 1-3 x 10 10 Mq based on a Salpeter mass function cut off at M = 0.1 M . 
The projected velocity dispersions are computed by a 2 = a 2 , cos 2 9 + cr 2 ,, sin 2 6, 
a 2 = cr 2 / sin 2 9 + cr 2 ; cos 2 9, and a z = o z i . The resultant velocity dispersions in the 
(x,y,z) coordinates are (a x ,a y ,a z ) = (110,82.5,66.3) kins" 1 . For the detailed 
computation of the optical depth and event rate based on the barred bulge model, 
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see Han (1994). The transverse velocity distributions for various models are listed 
in Table 1. 

There is disagreement over the three-dimensional distribution of stellar pop- 
ulations of the Milky Way. Therefore, it is useful to compare the effects on the 
optical depth, event rate, and reduced velocities using various models. We adopt 
three bulge models, isothermal, barred, and Kent (1992), and three disk models, 
Bahcall (1986), Kent (1991), and KP (1994). The density distributions of the mod- 
els are listed in Table 2. Since the mass-to-light ratio of the disk stellar population 
are not well determined, we normalize all disk models to the mass density of solar 
neighborhood po ~ 0.06 M Q pc~ 3 . Note that this value is somewhat higher than 
the observed stellar mass density near the Sun because the models are aimed at 
reproducing the star distribution ^ 1 scale height from the plane, not the local 
density. The density profile through Baade's Window at I — 1°, b — — 3°.9 is shown 
with respect to the distance from the Earth in Figure 1. 

Because various researchers have adopted different models, they find different 
optical depths and probabilities of lensing events from one another. To see the 
effects on the optical depths of different mass models, we compute the optical 
depths for various models and the results are listed and compared with values 
estimated by others in Table 3. The computed optical depths for various models 
are similar to each other and are not strongly model-dependent except for the 
barred bulge model: a barred bulge doubles the predicted bulge optical depth. 
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4. Lens Locations and Masses 



4.1. "Reduced" Velocity 

It would be very hard to untangle the contributions of bulge and disk Machos 
from the analysis of uo^ 1 alone. Griest et al. (1991) and KP (1994) point out that 
the mean duration of events uj^ 1 would be longer for lenses in the disk than for 
ones in the bulge. However, this occurs mainly because the event-rate distribution 
for disk events has a large tail toward the longer events. By contrast, the most 
common characteristic time of events is almost the same regardless of whether 
the lenses are bulge or disk Machos making it difficult to distinguish between them 
(Evans 1994). The value of u^ 1 for the most common events is ~ 7 (M/0.1 Mo) 1 / 2 
days, for both bulge and disk Machos. 

However, it would be possible to distinguish the bulge and disk Machos if 
their reduced velocities could be measured. The rates of lensing events by bulge 
and disk stars are computed from the equations in § 2 and are presented with 
respect to v y and v z as contour maps in Figure 2. We use the Kent and Bahcall 
models as representative models for the bulge and disk density distributions. All 
the computations hereafter are carried out with density and velocity distributions 
toward Baade's Window unless another window is specified. The isofrequency 
contours for a rotating bulge model are drawn with solid lines, while the dotted 
lines represent contours for a nonrotating bulge model. The rotation of the bulge 
does not dramatically affect the distribution of T. The event rates as functions of 
v y and R e are shown in Figure 3. Note that very few events generated by disk stars 
have v y < 0, whereas those generated by bulge stars are about equally likely to 
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have negative as positive v y . Therefore, the reduced velocity v could be used for 
partial separation of bulge from disk Machos. Recall, however, that measurement 
of v generally requires an additional satellite. 

Fortunately, the events generated by disk and bulge Machos can also be dis- 
tinguished by the reduced speed which generally requires only one satellite. In 
Figure 4, the rates of events with respect to v and R e are shown for bulge and disk 
Machos. Note that the plotted values of R e are for an assumed mass M — 0.1 M . 
For other masses R e oc M 1 / 2 . The rates with respect to v only are shown in Figure 
5 for various models listed in Table 1. For the spherical bulge model (and for either 
constant or linear disk velocity dispersion model) the majority of the events lensed 
by the disk Machos would occur in a region v < 500 km s _1 while bulge Macho 
lensing events would be found in a wide range of v (0-3000 km s^ 1 ). This means 
that a significant fraction of the bulge Machos could be separated from the disk 
Machos: events with v < 150 km s _1 are likely to be due to disk Machos whereas 
bulge Machos would be responsible for the events with v > 600 km s^ 1 . However, 
especially if the bulge has a bar-like structure it would be hard to distinguish disk 
from bulge Machos for 150 ^ v ^ 600 km s _1 , i.e. a substantial fraction of events 
[Fig. 5 (b)]. This predicament could be largely resolved by observing stars through 
other windows toward galactic bulge. The event rate as a function of v through the 
window at I = — 8°.66, b = — 5°.98 (Blanco & Terndrup 1989) is shown in Figure 5 
(d). 
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4.2. Mass Spectrum Of Machos 

The reduced velocity v determined from the parallax measurement of lensing 
events could tell one the physical parameters of the Machos which is very difficult 
to obtain from present observations. 

First, the approximate locations of disk Machos could be determined from the 
reduced velocity measurements. Assume for the moment that we could discriminate 
between disk and bulge Machos from the determined v. Figure 6 shows the event 
rates with respect to D \ for different values of v: the curves represent T(D {) for 
v = 50, 100, 150, and 300 km s _1 from the left to the right, respectively. Here, 
the event rates are computed with the linear disk velocity dispersion model. For 
each value of v, the most probable events occur at different values of D \ enabling 
one to determine the approximate locations of the disk Machos. The assumption 
that v discriminates between bulge and disk Machos is valid over different ranges 
of v for different windows. For example, for Baade's Window separation can be 
made for v ^ 150 km s _1 , while for the / = — 8°7, b = —6° window, separation 
requires v ^ 250 km s _1 . See Fig. 5. 

The next application of the reduced velocity measurement is that it could 
provide a strong constraint on the mass spectrum of Machos. Once again, let us 
assume that disk Machos could be distinguished from bulge Machos. The only 
constraint one has up to now for the estimate of Macho masses is the charcteristic 
time u' 1 . For demonstration, say that a particular event has cu^ 1 = 10 days, which 
is the most common value measured by the MACHO group. Figure 7 (a) shows the 
probability distribution of masses of a disk Macho constrained by a single piece of 
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information, oo^ 1 = 10 days, while the probability distribution constrained by the 
additional information, v, is shown in Figure 7 (b). The typical standard deviation 
of the distribution with two constraints is <J\ og M ~ 0.2, whereas o"i og M ~ 0.6 
just using the information of u^ 1 . The additional information of v significantly 
reduces the uncertainty of the mass estimate, and thus strongly constrains the 
mass spectrum of Machos. 

Unfortunately, it is more difficult to measure the mass spectrum of bulge Ma- 
chos. The reduced Einstein ring radius and velocity can be approximated as 



The ratio M/v can be determined from the measured value of R e and v, but one 
needs to know about v to obtain M. This can be done only by careful modeling 
of the velocity distribution of bulge Machos. 

4.3. Proper Motions 

For the lensing events in which a lens passes near the face of a star, the light 
curve deviates from its characteristic 3 parameter curve [eq. (1.1)] due to the 
differential magnification of the lensed star. By using this principle, one could 
measure the proper motion of a lens (Gould 1994a, Nemiroff & Wickramasinghe 
1994), which is defined by Q = v/D i. In Figure 8, we plot the rate of lensing events 
as a function of f2. The lensing events both by the disk and bulge Machos would 
occur most probably around Q ~ 1.1 Qq, where Qq = vq/Rq. Thus, it is unlikely 
that one could distinguish bulge Machos from disk Machos by measuring of VL. In 




(4GM/c 2 )(R 2 /D ls ) and v ~ v(R /D ls ). Then, 




(bulge Machos). 



(4.2.1) 
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addition, this measurement can be carried out for only ~ 2 (M/0.1 M )- 1 /2% f 



events. 



There would also be a detectable spectral line shift due to the rotation of a 
lensed star. By combining the line shift with to" 1 (from the overall light curve) 
and with the stellar rotation speed (obtained from the spectrum), the lens proper 
motion could be measured (Maoz & Gould 1994). Almost all the bulge stars 
observed by OGLE and MACHO have color V — I ^ 0.83 which are approximately 
K0 for main sequence stars and G5 for giants. Herbig & Spalding (1953, 1955) 
have shown observationally that there is a cutoff in the distribution of rotation 
velocities for evolved stars, which seem to be present near GO in all luminosity 
classes IV, III, II, and lb. Gray (1991) also provides a distribution of v ro t ■ In his 
distribution the drop in rotation of the giants takes place between GOIII and G3III. 
Since nearly all lensing events occur later than the cutoff, it would be difficult to 
detect the line shift due to the slow rotation of lensed stars. Moreover, even if 
the shift were detected, one must know the rotation speed in order to measure the 
proper motion. The rotation speed is difficult to measure when it is of order or 
less than the turbulent velocities of the stars. 

It is important to note, however, that one can uniquely determine v, M, and D i 
of Machos if both the proper motions and parallaxes are measured. The distance 
to a Macho is 



Since D os ~ 8 kpc, one can determine D \ from measured values of v and Vt. 
Once D i is known, the mass and transverse speed of a Macho can be uniquely 
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(4.3.1) 



determined. 



5. Conclusions 

The results of the computations of the optical depths and rates of lensing events 
with respect to various potentially measurable physical parameters are summarized 
as follows: 

1. The optical depths and rates of disk lensing events do not for the most part 
depend strongly on the models that have been proposed, based on star counts 
and stellar dynamics. However, both the optical depths and event rates would be 
considerably higher if the bulge is barred. 

2. The events lensed by bulge Machos could be distinguished from the ones lensed 
by disk Machos by measuring the reduced speed of Machos, which can be done 
with a single satellite. 

3. The approximate locations and masses of a significant fraction of Machos could 
be determined from the reduced speed of the Machos. Therefore, v would provide 
an important clue to the physical distribution and mass spectrum of Machos. 

4. It would be difficult to measure the proper motions of Machos seen toward the 
bulge either spectroscopically or photometrically. Even when lens proper motions 
can be determined, they do not discriminate well between bulge and disk Machos. 
However, if both the proper motion and parallax of an event are measured, then 
the mass, distance, and speed of the Macho are uniquely determined. 
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